Suppression of ablation in femtosecond double pulse experiments 
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We report the physical reasons of a curious decrease in the crater depth observed for long delays 
in experiments with femtosecond double pulses. Detailed hydrodynamic modeling demonstrates 
that the ablation mechanism is dumped when the delay between the pulses exceeds the electron-ion 
relaxation time. In this case, the interaction of the second laser pulse with the expanding target 
material leads to the formation of the second shock wave suppressing the rarefaction wave created 
^— «j ■ by the first pulse. The evidence of this effect follows from the pressure and density profiles obtained 

' at different delays after the first laser pulse. 
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During last decade, femtosecond laser systems have found numerous applications in different areas, such as laser 
material treatment, nano-optics, and surface analysis One of the particularly promising applications of ultrashort 
' (— I '! laser pulses is laser-induced breakdown spectroscopy (LIBS) 0], or a remote technique that can be used to instantly 
O i' analyze different materials. The advantages of ultrashort laser pulses were demonstrated in several LIBS experiments 
JL [ resulting in a considerably improved spectral resolution. To further increase the accuracy of the analysis, several 
' • special configurations were proposed enhancing the intensity of the plume emission. One of the promising schemes 
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is based on double pulse (DP) laser set-ups [JlBl- Iii addition, experiments with several delayed femtosecond pulses 
^ ' were shown to be advantageous in laser machining Q , thin film modification 0] and in fast ion production ^ . 

• [ Performing DP femtosecond experiments for different metals, several authors 0, 0, 0] surprisingly noticed a mono- 
Q ■ tonic decrease in the resulting crater depth with the increase in the delay rngiay between the laser pulses. This effect 
was observed both under vacuum 0] and in air Q . In the experiments [5| with copper targets, for instance, the 
J>Lj' laser fluence of each pulse was set to be Fgingie = 2 J/cm^. When the delay between the pulses was much shorter 
,Xh I than the electron-ion relaxation time Tei in the target material (rdciay *C r^i, where r^i « 10 ps for copper [13] ) the 
^ crater depth was the same as in the case of a single pulse (SP) with the laser fluence 2 x i^singie = 4 J/cm^. For the 

delays approaching the electron-ion relaxation time (rdoiay ~ ^ei), the crater depth monotonically decreased. Finally, 
fSJ ■ for longer delays {r^cUy 3> r^i) the crater depth was found to be even smaller than that obtained with a SP at the 
^ ' fluence i^singio, see Fig. [TJ Similar behavior was also obtained for aluminum 0], gold 0| and nickel ^ thus indicating 
■ the common physical effect. The dashed (blue) curve in Fig. [T]is just a prediction for long delays, which have not 
been yet achieved in the DP experiments. This case will be discussed at the end of the letter. 

Recently, several explanations of the unusual dependency of the crater depth on the delay were proposed. In 
_^ _ particular, the temperature dependence of the heat conductivity claimed to explain the observed effect in [5|. Target 
. heating can in fact lead to the energy accumulation changing target absorption properties. In addition, the energy 
I absorption in the laser-generated plasma plume was furthermore evoked in [3|. The latter explanation makes sense if 
. the time delay is sufficiently long, so that the second pulse interacts with the laser-created expanding plasma plume. 
Recently [§],the ablation reduction effect was also attributed to the presence of a liquid layer [nf . In fact, laser 
ablation of metals is commonly accompanied by melting. The supercritical concentration of the conduction electrons 
in liquid phase makes, however, the melted layer non-transparent for laser light. Therefore, a better understanding of 
the experimental findings 0, |^, 0| requires a more complex self-consistent analysis of an interplay of different physical 
effects to highlight the basic mechanism of suppression of ablation. A detailed hydrodynamic modeling can provide 
a comprehensive physical scenario of the interaction dynamics. In the present letter, we report the results of such 
numerical calculations based on a two-temperature hydrodynamic model. The model accounts for such processes as 
laser light absorption, plume expansion, thermal conductivity, two-temperature effects, phase transitions, material 
decomposition, etc. In particular, the dynamics of the DP ablation is investigated in details and is shown to play a 
crucial role in the considered phenomenon. 

The hydrodynamic model [l2,ll3] was developed and previously used for the investigation of a SP ablation of metals. 
The main advantages of the model are in the treatment of mctastablc phases and in the used multiphase wide-range 
equation of state [l3| . The validity of the model was confirmed by the comparison with the experimentally-obtained 
ablation depth for several metals [IJ]. In the case of a DP, the model accounts for the laser energy absorption by the 
ablation plume generated by the first pulse. To treat the absorption process, we use the Helmholtz wave equation with 
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the quasi-classic solution described elsewhere 15[. The absorbed laser power density (per unit volume) is expressed 
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in this case as Ql = IokQlm.{e} \E/Eo^ . Here /q is the peak laser intensity, fcg = ^l/c is the laser wave vector 
with the laser frequency wl and the light speed in vacuum c, e is the complex permittivity of media, x) is 
the solution to the Helmholtz equation and Eq = ^/SttIqJc is the maximum of laser electric field in vacuum. The 
wide-range model of permittivity [lit describes Drude-like (intraband) effects, band-to-band transitions and the hot 
plasma limit. The electron thermal conductivity in metals can be measured experimentally only for low and near- 
room temperatures [l7l|. T o extend the model jl8| used in this work over the range of high temperatures we apply 
the first-principle data 1^, 23). The model for electron- ion coupling (electron-phonon for low temperatures) is similar 
to one used in [lit ^^'^ is also adjusted to satisfy the VASP ( Vienna Ab-initio Simulation Package) first-principle 
calculations (2l|. 

Calculations arc performed for copper with two s-polarized Gaussian laser pulses of 100 fs width, 800 nm wavelength, 
and with first and second pulse intensity maximum at t = and t = Tdeiay, respectively. Fig. [5] shows the dynamics of 
both shock and rarefaction waves in the copper target during and after the propagation of the first pulse. At < = ps, 
the ion pressure increases due to the collisions with hot electrons. After the electron-ion collisional relaxation, the 
temperatures of the electron and ion sub-systems are the same, but the pressure of ions is much larger (up to 35 GPa 
at 10 ps in Fig. ^jp). At the moment of 20 ps after the peak of the first pulse, one can see the appearance of the 
negative pressure zone (with pressure about —2.5 GPa at a; « 35 nm). Later, at moments of 50 ps and 75 ps after 
the maximum of the first pulse, the negative pressure value is about —3.6 GPa located at x « 100 and x w 175 nm 
depth, respectively. Such pressure values were shown to result in the mechanical spallation in the liquid phase, but 
are insufficient for the spallation in the solid phase. The arrows in Fig. [2^ show the position of the melting front 
x^cit(t) for the corresponding moments. It is seen that the melting front propagation speed drops as the laser energy 
dissipates. For x > Xmdt the substance is in solid phase and the spallation strength for a solid material is close to its 
theoretical limit [S^l . 

Our previous calculations 12:ll3| demonstrated that only the melted region can be ablated. Therefore, if the second 
pulse arrives during the rarefaction wave propagation (10 ps < i < 50 ps in Fig. |2^, b) through the liquid layer, this 
pulse can reheat the nascent ablation plume. As a result, a high-pressure region is generated in the vicinity of the 
initial target surface. This hot region then produces a shock wave and thus diminishes the action of the rarefaction 
wave. 

Fig. [3] illustrates the location of the laser energy absorption for both the first and the second pulses. The ratio 
of the free electron concentration rie to the critical one Ua = w|_me/(47re^) is presented in Fig. |3K. When the ratio 
rie/nci- > 1, one can observe the common skin effect for the first pulse. For the second pulse, the electron concentration 
gradient is smooth and the absorption region is shifted to the left of the initial free surface position a; = nm. The 
spike in the electron concentration in Fig. [3^ at 50 ps correlates with the broken off chunk seen in the density profile 
in Fig. for the same moment. Thus, the absorption of the first pulse takes place from to 50 nm in the skin 
layer (solid (red) curve), whereas the second pulse is absorbed from —200 to nm (dash-doted (blue) curve) with the 
maximum at about —100 nm (rie ~ 0.03ncr)- 

To clearly illustrate the suppression of the rarefaction wave in the DP experiments. Fig. 2] presents the pressure 
evolutions calculated at the depth of 100 and 200 nm under the initial surface position. For the delay rjoiay = ps 
(Fig. 2^), the amplitude of the rarefaction wave is about —3.5 GPa both in the calculations with SP and DP. The 
melted layer can, therefore, be ablated forming a crater. At the delay of 5 ps, the rarefaction wave is formed later, 
and the amplitude of the wave decreases (~ —3 GPa). For longer delays, the minimal pressure appears also later 
(at 200 nm depth and for 10 (Fig. and 25 ps delay the moments are 100 and 125 ps, respectively), whereas the 
amplitude of the negative pressure decreases (at 200 nm depth for 10 and 25 ps delay the pressure values are —2.8 
and —2.3 GPa, respectively). The observed attenuation of the pressure amplitude results in the decrease of the crater 
depth. This effect is observed despite the fact that in our simulations the the melted layer is thicker for longer delays 
between pulses. Finally, for the delays of 50 and 100 ps (Fig. Hj;, d), the first pulse leads to the spallation of the 
melted layer, whereas the intensity of the second rarefaction waves (at about 150 and 200 ps, respectively) is about 
— 1.8 GPa. In this case, the ablated depth is close to that obtained for the SP with 2 J/cm^. 

The performed simulations give us an opportunity to observe the dynamics of the ablation process. We suppose 
here that the crater depth is formed due to the material removing after the laser irradiation. We can calculate 
the ablation depth Aabi integrating the mass flux through the plane x = (initial target surface) by using the 
equation Aabi(i) = Pa^ jQ{pu)\x=odt' , where po is the initial material density and u is the material speed. As one 
can see in Fig. [SI the removed mass drops while the delay between pulses increases. Already for the 10 ps delay the 
crater depth for a DP case is equal to that obtained for a SP (short dotted (green) curve and solid (black) curve in 
Fig. \5\ respectively). For longer delays (50 and 100 ps) the crater depth is even smaller than in the case of a SP. 
The absorption of the second pulse in the nascent ablation plume results in the reheating of ablated material and 
acceleration of outward part of the plume and deceleration of the inner part of the plume. One can see in Fig. [5] 
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(ideiay = 100 ps) that the back flux dominates from 170 to 280 ps and results in temporal decrease in the crater depth 
for 100 ps delay. The corresponding experimental points arc presented to prove the validity of the model. 

We can also suggest the crater formation mechanism in the case rdciay — > oo (Fig. [U dashed (blue) curve) . During 
the plume expansion, the electron concentration drops and the plasma cloud becomes transparent for the second 
pulse. This effect can be achieved for itransp ~ R/u with the radius of the plume R and the rarefaction velocity u. 
Simple estimations yield itransp on the order of hundreds of nanosecond. 

In summary, we have revealed the main stages of the DP ablation process. The first pulse is absorbed by the 
conduction band electrons in metal leading to a rapid rise in their temperature and pressure. During the electron-ion 
relaxation, the temperature of ions reaches the melting point, so that homogeneous melting occurs. At this time 
(t ^ Tf^i), the pressure of ions rapidly grows resulting in the formation of a shock wave propagating into the target. 
At the same time, a rarefaction wave is formed near the free surface. This rarefaction wave propagates through the 
melted surface layer and gives rise to the mechanical fragmentation and ablation of the liquid target material. In the 
DP experiments, the delay between laser pulses plays a fundamental role. When the delay is much shorter than the 
relaxation time, only one shock wave and one rarefaction wave appear. In this case, the ablation crater is formed by 
both pulses simultaneously as in the case of a SP of the same energy. When the delay is on the order of r^i , the second 
pulse reduces the intensity of the rarefaction wave and the depth of the ablation crater decreases. Finally for long 
delays Tdeiay ^ Tei the ablation crater is formed by the first pulse only, whereas the second pulse reheats the ablated 
material and does not cause any additional fragmentation of the target. We conclude that the demonstrated effect of 
the suppression of the rarefaction wave by the second laser pulse explains clearly the results of DP laser experiments. 
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FIG. 1: (Color online) Empty (red) circles — DP experiment 0] with Cu ablation by pulses with -Fsingio 
curve — interpolation of the experiment. Dashed (blue) curve — our theoretical assumption for Tdday oo- Time Tt 



2 J/cm^; solid (red) 
is the 



moment of the ablated substance transparency for the second pulse due to rarefaction. Horizontal dash-and-dot (black) lines 
show depths from SP and two pulses separated by rdeiay — >■ oo. 
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FIG. 2: (Color online) Calculation results for Cu target ablation in vacuum by a SP with fluence 2 J/cm^. Temperature — (a), 
pressure — (b) and material density — (c) of the ions are shown for different moments. The initial free surface of the target 
is located at a = nm, the laser pulse goes from the left. Arrows show the melting front position a::mcit for different times 
(1—0 ps, 2—10 ps, 3—20 ps, 4—50 ps, 5—75 ps). 
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FIG. 3: (Color online) DP with rdeiay = 50 ps for Cu target. The ratio of the electron concentration to the critical one — (a), 
relative absorption efficiency — (b). The first pulse maximum is at t = ps, the second one is at t = 50 ps. Initial target 
position is at a; = nm. 
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FIG. 4: (Color online) Evolution of the ion pressure for 100 and 200 nm depth in Cu target. Results presented for different 
delays between femtosecond pulses 2x2 J/cm^: solid (red) curve — 100 nm depth, long dashed (blue) curve — 200 nm depth. 
SP pressure profiles for 2 J/cm^ are also presented: dash-and-dot (green) curve — 100 nm depth, short dashed (magenta) 
curve — 200 nm depth. 
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FIG. 5: (Color online) Ablation of Cu target. Simulation of crater depth Aabi growth in time for different DP delays. Crated 
depth from experiment for different delays is also presented: filled (red) circle — ps, filled (magenta) square — 5 ps, filled 
(green) triangle — 10 ps, empty (cyan) circle — 50 ps, empty (blue) triangle — 100 ps, empty (black) square — SP. 



